subroutine r2abij_t2efij_iefab_oc_dg(r2, t2, &
              iabcd, icore, fact)
!
use mod_iop
use mod_size
use mod_orbit
use mod_ioff 
implicit none
!
real*8, intent(in) :: fact
real*4, intent(in) :: iabcd(*)
real*4, intent(in),device :: r2(*)
!
real*4, intent(inout),device :: t2(*), icore(*)
!
integer :: i, j, a, b, nsize, i0, i1
real*8, external :: dnrm2
!this is to do r2 = r2 + 0.5*t2(efij)*iefab
i0 = 1
i1 = i0 + nvvoo*isd
!call printdnrm2('iabcd', iabcd, 6)
! ABAB  
call VmnpqVmnrs_to_Vpqrs_large_g2(t2(ioi3(2)), icore(i1), 1, 1, vrta, vrtb, 1, &
               popa, popb, 1, vrta, vrtb, 1,  icore, 1, 1.d0, & 
               0.d0,'ABCDAB__',iabcd)
call itranspso(icore, icore(i1), popa, popb, 1, vrta, vrtb, 1, 1) 
call axpyso(nioi32, fact*2.d0, icore(i1), 1, r2(ioi3(2)), 1)
!
return
end
